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ABSTRACT 

The recent detection of a 7-ray flux from the direction of the Galactic center by 
EGRET on the Compton GRO raises the question of whether this is a point source 
(possibly coincident with the massive black hole candidate Sgr A*) or a diffuse emitter. 
Using the latest experimental particle physics data and theoretical models, we examine 
in detail the 7-ray spectrum produced by synchrotron, inverse Compton scattering 
and mesonic decay resulting from the interaction of relativistic protons with hydrogen 
accreting onto a point-like object. Such a distribution of high-energy baryons may be 
expected to form within an accretion shock as the inflowing gas becomes supersonic. 
This scenario is motivated by hydrodynamic studies of Bondi-Hoyle accretion onto 
Sgr A*, which indicate that many of its radiative characteristics may ultimately be 
associated with energy liberated as this plasma descends down into the deep potential 
well. Earlier attempts at analyzing this process concluded that the EGRET data are 
inconsistent with a massive point-like object. Here, we demonstrate that a more careful 
treatment of the physics of p-p scattering suggests that a ~ 10 6 M Q black hole may be 
contributing to this high-energy emission. 
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1. Introduction 

X-ray and 7-ray emission have been detected 
from the Galactic center flWatson et al. 1981 



Skinner et al, 1987] ; [Predehl Sz Triimper 199*4 ; 



Chnrazov et al. 1994] ) . The implications for the 
radio point source Sgr A* are rather interest- 
ing, since the X-ray luminosity is not as large as 
what is expected based on X-ray observations of 
smaller black hole candidates. Recently, EGRET 
on board the Compton GRO has identified a 
central (< 1°) ~ 30 MeV - 10 GeV continuum 



source with luminosity ~ 5 x 10 36 ergs s 1 ( [Mat- 



box 1996| ; |Mattox et al. 1996| ; [Merck et al. 1996[) 
This EGRET 7-ray source (2EGJ1746-2852), ap- 
pears to be positioned at / as 0.2°, but the exact 
center of the Galaxy, or even a negative longi- 
tude, cannot be ruled out completely. Its spec- 
trum can be fit by a hard power-law of spectral 
index a = -1.74 ± 0.09 {S = S E a ), with a 
cutoff between 4 — 10 GeV. At lower energies, 
the COMPTEL data provide useful upper limits 
( (Strong 1996|) . 

These 7-rays may originate either (1) close to 
the massive black hole, possibly from relativistic 
particles accelerated by a shock in the accret- 
ing plasma ( Mastichiadis fc Ozernoy 1994] ), or 
(2) in more extended features where relativistic 
particles are known to be present (Pohl 1997). 
In the former (see also Mahadevan et al. 1997 . 
who considered a thermal distribution of hot pro- 
tons), the 7-rays may result from the decay of 
7r's produced via p-p collisions of ambient pro- 
tons with the shock-accelerated relativistic pro- 
tons. This study concluded that the presence of 
a Mh ~ 10 6 Mq black hole was inconsistent with 
the EGRET data. However, the earlier calcula- 
tions suffered from over-simplification and an in- 
complete treatment of the physics of p-p scatter- 
ing. For example, although the multiplicity of tt 
production (i.e., the number of 7r's produced per 
collision) is a strong function of energy, it was 
approximated with a constant value of 3 when 
in fact it can change by orders of magnitude. 



In addition, ignoring the role of cascading pro- 
tons is not a valid approximation when the en- 
ergy carried away by the leading protons can be 
over 90% of the incoming proton energy. Here we 
examine the hypothesis of a black hole origin for 
the 7-rays employing the most current data for 
the energy-dependent cross-sections, inelasticity, 
and tt multiplicity, together with a self-consistent 
treatment of the particle cascade. 

2. The Physical Picture 

Sgr A* is a nonthermal radio source at the cen- 
ter of the Galaxy (e.g., Menten et al. 1997]). The 
large mass (~ 1.8 x 10 Mq; |Haller et al. 1996] ; 
Eckart & Genzel 1997) enclosed within ~ 0.1 pc 
indicates that this object may be a massive black 
hole. The nearby IRS 16 cluster of hot giant 
stars produces a Galactic center wind with ve- 



locity v 



gw 



500 — 700 km s 1 and mass-loss 



rate M gw 3 - 4 x 10~ 3 M yr' 1 . A por- 
tion of this wind is captured by Sgr A* and 
accretes inward. The accretion rate ( ^ 10 22 g 
s _1 ) resulting from this Bondi-Hoyle process is 
well below the Eddington value, and so the gas 
attains a free-fall velocity ( [Mclia 1994 ; Ruffert 
fc Melia 1994] ), which eventually becomes super- 
sonic, and a shock forms at r s h ~ 40 — 120 r g 
flBabul et al. 1989| ), where r g = 2GM h /c 2 is the 
Schwarzschild radius. A fraction of the particles 
may be accelerated to very high energy by the 
shock. However, the greater synchrotron and in- 
verse Compton efficiency of the e~'s compared 
to that of the p's limits the maximum attain- 
able Lorentz factor of the former by several or- 
ders of magnitude compared to the latter, and so 
the (accelerated) e~ contribution to the radia- 
tion field is negligible. The relativistic p's are 
injected through the shock region with a rate 
Pp(Ep) = p E~ x cm -3 s -1 GeV -1 . In steady 
state, this leads to a power-law distribution with 
index z ~ 2.0 - 2.4 ( pones fc Ellison 1991|) . In 
our case, z is determined in part by the p cool- 
ing processes and the particle cascade, and as we 
shall The normalization p Q is related 



2 



to the efficiency 77 of the shock by 



J PoE p X E P dEp — r]L gr , 



r]GM h M 
r s h 



(1) 



These relativistic particles interact with the 
ambient particles and the magnetic field B, pro- 
ducing photons via synchrotron, inverse Comp- 
ton scatterings and the decay of mesons created 
during p-p collisions. Because the injected p's 
can be ultrarelativistic, leading order nucleons 
produced in the scattering events also contribute 
to the spectrum via multiple collisions in an en- 
suing cascade. The main products in these col- 
lisions are 7r's, which then decay either to pho- 
tons (tt° — ► 77) or leptons (jt — > fx ,Up, with 
fj,^ -> e ± u e Vp). 

Following Melia (1994), one can see that on 
average the ambient p number density is 



N 



n. 



47rcr| 



3/2 



(2) 



where N ~ M/m p is the p number accretion rate 
in terms of the mass m„. If in addition the mag- 
netic field B is in approximate equipartition with 
the kinetic energy density, then 




(3) 



The frequency v max at which the gas becomes 
transparent lies in the Rayleigh-Jeans portion of 
the spectrum. This radiation is assumed to be in 
thermal equilibrium with the e~'s. 

3. Properties of the Particle Cascade 

The relativistic p 7 s undergo a series of inter- 
actions including pN — ► pN M v M WJ y, where N 
is either a p or a neutron n, M n represents the 
energy-dependent multiplicity of it's, and M N ft 
is the multiplicity of nucleon/anti-nucleon pairs 
(both increasing functions of energy). Since 
M NR /M n < 1CT 3 at low energy and even smaller 



at higher energies ( pline 1988 ), we here ignore 
the anti-nucleon production events. The charge 
exchange interaction (p — > n) occurs roughly 
25% of the time (e.g., |Begelman et al. 1990| ). 
The other possible interactions are ^7 — ► pvr°7, 
P7 — ► n-7r + 7, p7 — > e + e~p and pe — ► eNM n . 

The high-energy cutoff for the injected p dis- 
tribution is set by determining the Lorentz fac- 
tor "fp^max above which the combined energy loss 
rate due to synchrotron emission, inverse Comp- 
ton scattering and hadronic collisions exceeds the 
rate of energy gain due to shock acceleration. 
This transition energy depends on the functional 
form of the inelasticity and the fraction of power 
transferred to the 7r's during the p-p collisions. 

Using the M n measured at several center-of- 



momentum (CM) energies ( Alper et al. 1975 ; 
Abe et al. 1988; Albajar et al. 1990), one can 



determine the ir injection rate from the p dis- 
tribution and the physical characteristics of the 
ambient medium. The particle cascade contin- 
ues with the emission of 7-rays and leptonic de- 
cays. The e~'s and e + 's produced in this fash- 
ion constitute an energetic population and one 
must assess their contribution to the spectrum 
via synchrotron emission and inverse Compton 
scattering. 

For the conditions in Sgr A*, the p-p collisions 
dominate over all other it production modes. The 
relevant energy (E p ) range is bounded below by 
the 7r production threshold and above by 7 Pjma x- 
Using logarithmic bins, and assuming time inde- 
pendence, we first calculate the steady state p 
distribution p p (E p ) using the diffusion loss equa- 
tion 



Pp(E p 



k„ E, 



p p dE 



.892 J R P p{Ep)AdE p + pp{Ep) 

I [n p a pp (Ep)c + 2E p k p ] , (4) 



2 d Pp 



where A 



and 



5{E' p -E p>1 ) + 5{E' p -E p>2 ) 

Rpp(Ep) = n p a pp (E p )cpp(E p ) cm~ 3 s~ 1 GeV~ 1 is 
the rate of p-p collisions at energy E p , and k p = 



3 



(AaTTnl/3c 3 mp)(UB + U^) is the sum of constants 
appearing in the power P syn c+Pcom P ton « k p E%. 
The first term represents the influx of p's hav- 
ing energy E p as secondaries in p-p collisions 
between relativistic p's with energies E p i and 
S Pi 2 and the ambient p's (which are effectively 
at rest). The relationship between E p and E p \ 
and E p> 2 is unique, as expressed by the delta- 
functions, and is determined by special relativity 
and the inelasticity K pp . We use the assump- 
tion that on average the two leading p J s created 
travel either parallel or opposite to the boost to 
simplify our calculations. In the p-p CM frame, 
we use K pp = 1.35 s" - 12 for > 62 GeV, and 
K pp = 0.5 for ^/s < 62 GeV, where the higher en- 
ergy slope is from Alner et al. (1986), normalized 
to match the approximately constant low energy 



value ([Fowler et al. 1984| ). 

The cross section a pp is taken as a function 
of energy from the most current published data 



( [Barnett et al. 1996 ). Unfortunately, the high- 
est energy achieved in modern colliders is orders 
of magnitude below the values attained in our 
system. However, the data for y^s ^ 100 GeV 
have a log-linear form which makes it possible 
to extrapolate up to much higher E p . For the 
entire range, this is within the Froissart upper 
bound, which states that at extremely high en- 
ergy, (T P p oo oc (Ins) 2 . 

The steady state relativistic p distribution re- 
sulting from this procedure can then be used to 
calculate the synchrotron and inverse Compton 
scattering spectra, following Rybicki & Lightman 
(1979), and the rate R pp of p-p collisions. For 
each of these collisions, a multiplicity M v of 7r's 
is produced, with a ratio of charged to neutral 
particles of roughly 2:1 (exactly when there are 
two leading p's, otherwise there will be a surplus 
7r + to conserve charge). These have a distribu- 
tion in transverse (to the beam in experiments, in 
our case to the direction of the boost back to the 
lab frame) momentum dN n /dp±, which is mea- 
sured as a function of -y/s at collider experiments. 
In order to find the energy of the 7r's in the CM 



frame, we also need the parallel component of the 
momentum, p» , which we extract from the ir dis- 
tribution as a function of the rapidity, y. In the 
CM frame, y = (1/2) ln[(£;+ P ||)/(£* -p,,)], and 
y ~ — ln[tan(0/2)] for relativistic energies, where 
cos# = pii/ | p \. At lower energy (y/s ^ 200 
GeV), dN n /dy is Gaussian in shape, the top of 
which widens gradually into a plateau with in- 
creasing energy. The width and the height of 
this plateau can be fit to functional forms in 
yfs. Details of this process are discussed in an 
upcoming paper. After binning in p±, the end 
product is a distribution dN 7T /dp±dy of it's in 
the CM frame, for each p-p collision at any par- 
ticular CM energy. Given y and p», and the 
fact that E* = (m 2 + p\ + P±) 1/2 , the dis- 
tribution (cm _3 s _1 per p-p collision at E p ) of 
7r's at energy E* follows from a convolution of 
(dN n /dp±dy)dp±dy with R pp . 

Each of the photons produced in the ir° decay 
acquires a rest frame energy e 7 = (l/2)m 7r c 2 and 
is emitted with equal probability in any direction. 
The photon number density in the observer's (or 
lab) frame must therefore be 



p 7 (e 7 ) de~ 



2de, 



(5) 



where 7^ and (3 n are the tt Lorentz factor and di- 
mensionless velocity, respectively, in this frame. 
Finally, the photon emissivity is given by the 
multiple integral expression 



j 7 (e 7 ) = / p 7 (e 7 ) 



-R pp (Ep)dE p dp± dy 



dp±dy 

photons cm 3 s -1 MeV -1 . (6) 



The charged it's decay to leptons, which can 
themselves be a source of radiation from syn- 
chrotron and Compton processes. We follow a 
procedure for the e ± completely analogous to 
that developed above to find the steady state 
distribution: 



Pe{E e 



-2k £ E e p e (E e ) - K e E\ 



dp e {E e 
dE e 



(7) 
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This equation is easy to solve for the distribution 
p e (E e ), which we then use to calculate the syn- 
chrotron and inverse Compton spectra from the 
cascade e ± s. Here, n e = (4/3)<JTc/(m e c 2 ) 2 (UB + 

4. Results and Discussion 

When the shock is located at r s h = 40r s , we 
have B ~ 435 Gauss and n p ~ 1.6 x 10 9 cm~ 3 . 
The Rayleigh-Jeans tail cuts off at v max « 10 13 
Hz, with a temperature ~ 6.3 x 10 9 K. Thus, 
lp,max ~ 1-2 x 10 9 . We find that for the environ- 
ment of Sgr A*, unlike those of typical AGNs, 
the p-p collisions dominate over p-j by at least 
ten orders of magnitude over the entire energy 
range. This is because of the relatively low den- 
sity of ambient baryons, and the extreme dearth 
of photons due to the low value of v max . Because 
of this, we neglect the contribution from all p-7 
interactions. We can similarly neglect p-e colli- 
sions for this system, the cross-section of which 
is igyOpp- Since the population of relativistic 
e~'s is much smaller than that of ambient p's, the 
contribution is insignificant. The synchrotron 
cooling channel dominates at the highest energy, 
and is the other significant contributor to the 
spectrum. From Sikora et al. (1989), we find that 
any relativistic n's with Lorentz factor 7„ ^ 10 8 
will escape the system without interacting, so we 
consider as lost any n's produced in the cascade. 
The low photon density also means the region is 
extremely optically thin to high-energy photons, 
so we consider the spectrum observed at Earth to 
be that produced at the shock. By comparison, 
le,max ~ 6.5 x 10 5 for the e~'s, corresponding 
to an energy of 3.3 x 10 5 MeV. For a roughly 
equal injection rate of relativistic e~'s and p's, 
this energy content is significantly below that of 
the cascade e~'s, which therefore dominate the 
leptonic contribution to the photon spectrum. 

In all, five spectral components may be con- 
tributing to 2EGJ1746-2852: p synchrotron, p 
inverse Compton scattering, e synchrotron, 
inverse Compton scattering, and 7r° decay. For a 



reasonable efficiency (i.e., r\ ^ 10%), the p syn- 
chrotron spectrum dominates over that of the 
photons from ir° decay as long as the p injection 
index x ^ 2.2. In Figure 1, we show these compo- 
nents for the case when r s /j = 40r 9 , and x = 2.0 
with an efficiency of 1%. The p synchrotron 
seems to fit the data reasonably well, but clearly 
misses the apparent low energy turnover in the 
EGRET data, and the upper limits for the high- 
est energy COMPTEL points. This could be 
due to the simplified geometry we have adopted 
in this paper, but this is unlikely since the syn- 
chrotron spectrum depends primarily on the par- 
ticle physics. It is also evident in Figure 1 that 
Compton scattering is not important for this 
source, and that the cascade e ± 's are relatively 
ineffective. The photons produced by tt decay 
are also insignificant. Placing the shock at 120r 9 
instead of 40r 9 (Fig. 2) increases the emission 
area while decreasing B, and the p synchrotron 
spectrum misses more data at the low end of the 
EGRET spectrum. The components still do 
not contribute. 

With the shock at 40r 5 , the tt° decay spectral 
component begins to dominate over synchrotron 
when x ^ 2.2. In Figure 3, we show the same five 
components for the case x = 2.4, which leads 
to z = 2.46. Here, r\ ~ 9%. The shape of the 
7r°-induced 7-ray spectrum can be understood as 
follows. The center of the curve is set by the en- 
ergy (e 7 = 67.5 MeV) of the decay photons in 
the 7T° rest frame, and the width is determined 
by Doppler broadening. The slope of the sides, 
and hence the index of the EGRET spectrum, 
is due to the falloff in the number of decaying 
7r's at higher energy. Each tt decay produces a 
flat photon spectrum (see Eq. |5|) whose width 
increases with E n . So the cumulative effect of 
all the decays is greatest near e 7 = 67.5 MeV, 
where all 7r's contribute. The relative contribu- 
tion to the spectrum at lower or higher e 7 then 
depends on the overall tt distribution, which in 
turn is a function of both M n and p p (E p ). The 
flattened top, and hence the low-energy turnover 



5 



in the tt decay spectrum, is due to the cutoff in 
tt production near the threshold. The tt decay 
spectrum cannot be translated laterally, so a si- 
multaneous match of both this turnover and the 
spectral slope is significant. 
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Fig. 1. — Five spectral components (as labeled) 
resulting from the p-p cascade. The shock is here 
located at 40r 9 , and the injected proton spectral 
index is x = 2.00, yielding a steady state proton 
index z = 2.06. The inferred efficiency in this 
case is 77 = 1%. The data points are from: (radio) 
Lo (1987), Zylka et al. (1993); (IR) Eckart et al. 
(1993), Stolovy (1996); (X-rays) Pavlinskii et al. 
(1992); (7-rays) Strong 1996, Mattox (1996). 

The required efficiency for producing the EGRET 
spectrum is different depending on whether syn- 
chrotron or tt decay emissivity dominates. There 
are three exit channels for the processed p en- 
ergy: (i) p synchrotron, (ii) tt° — ► 77, and (iii) 
tt^ — * fJ, Vn, with / u ± — ► e^UeU^. The latter two 
are coupled by a fixed ratio since each of the three 
types of 7r's receives an equal fraction of the en- 
ergy lost by the colliding p's. When x 2.2, (i) 
dominates and the conversion of p power to pho- 
tons is very efficient. When x £ 2.2, (ii) & (iii) 
dominate, but the fraction of p power going into 
photons rather than particle by-products is now 
smaller, so a higher 77 is needed. 
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Fig. 2. — Same as Fig. 1, except for a shock 
located at 120r 9 . 

A population of relativistic p's energized within 
an accretion shock near a super-massive black 
hole at the Galactic center may be contribut- 
ing to the ~ 30 MeV-10 GeV emission from this 
region. Depending on the value of the injected 
p index (i.e., x ^ 2.2 or x £ 2.2), this contribu- 
tion may come either from synchrotron or tt° de- 
cay. However, the synchrotron emissivity cannot 
account for the turnover in the EGRET spec- 
trum and some of the COMPTEL upper limits 
at e 7 ~ 100 MeV, whereas the 7r-induced photon 
distribution has a natural flattening there due to 
the threshold for tt production in p-p scatterings. 
We have not yet undertaken a detailed x 2 fitting 
to find the optimal parameters. However, it may 
be an indication of robustness in the model that a 
good fit was obtained without any fine-tuning. A 
more sophisticated fitting using the actual data 
will be undertaken in the future. Our methods 
and results differ significantly from earlier treat- 
ments. Our use of M n and a careful treatment 
of the particle cascade give a good fit to the 7- 
ray spectrum with reasonable parameters, such 
as z ~ 2.5, 77 ~ 0.09, and a black hole mass 
M h ~ 10 6 M (cf. Mastichiadis fc Ozernoy 1994j 
who concluded that z « 1.7, and M h < 10 6 M ). 



6 
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Fig. 3. — Same as Fig. 1, except that the spec- 
trum is here magnified to highlight the EGRET 
and COMPTEL energy range. In addition, the 
proton injection index is x = 2.4, yielding a 
steady state proton index z = 2.46. The inferred 
efficiency is here ij = 9%. 
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